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Complex networks characterized by global transport processes rely on the presence of directed 
paths from input to output nodes and edges, which organize in characteristic linked components. 
The analysis of such network-spanning structures in the framework of percolation theory, and in 
particular the key role of edge interfaces bridging the communication between core and periphery, 
allow us to shed light on the structural properties of real and theoretical flow networks, and to 
define criteria and quantities to characterize their efficiency at the interplay between structure and 
functionality. In particular, it is possible to assess that an optimal flow network should look like a 
"hairy ball" , so to minimize bottleneck effects and the sensitivity to failures. Moreover, the thorough 
analysis of two real networks, the Internet customer-provider set of relationships at the autonomous 
system level and the nervous system of the worm Caenorhabditis elegans -that have been shaped by 
very different dynamics and in very different time-scales-, reveals that whereas biological evolution 
has selected a structure close to the optimal layout, market competition does not necessarily tend 
toward the most customer efficient architecture. 



I. INTRODUCTION 

Despite profound differences, natural and artificial net- 
worked systems share striking similarities. Complex net- 
works science [TJ [21 13] has successfully rationalized sev- 
eral of the most ubiquitous features, such as the small 
world property or the presence of strong degree hetero- 
geneity, relating them to the existence of general orga- 
nizing principles. These self-organization laws also shape 
the observed large-scale connectivity layout of a special, 
yet common, class of networks describing transport pro- 
cesses, be it of matter, energy, or information. These 
networks are characterized by asymmetric interactions 
giving rise to local flows that collectively organize into a 
large-scale stream dominated by a processing core which 
transfers input into output: the universal bow-tie archi- 
tecture |1] that is intimately related to the functional 
activity of these systems. 

In general terms, most previous research exploring the 
relation between form and function in complex networks 
has been mainly focused on the analysis of topological 
features such as modular ordering revealing functional 
aspects with fewer exceptions treating directly func- 
tional aspects such as efficiency [B]. Specifically, trans- 
port has been studied as one of the main functions influ- 
enced by topology fT*, 'B^ and functional design principles 
of global flux distributions have been discussed for biolog- 
ical networks 0113 EI]- Despite these efforts, the "form 
follows function" assertion still remains to be fully un- 
derstood from a complex networks science perspective, a 
major difficulty in the fact that present network patterns 
are the result of non-stationary and adaptive evolution- 
ary histories that can greatly vary for different networks. 
However, general self-organization principles should not 
only govern structure but also their interplay with func- 
tional features. 

Our purpose of inferring information about function 
and evolution from a precise knowledge of the topologi- 



cal makeup requires the understanding of how flow net- 
works organize to develop functionality. In this respect, 
percolation theory on complex networks |12j provides a 
valuable framework to discuss their connectedness and to 
identify the components that are key to a complete de- 
scription of their global connectivity layouts conforming 
the percolation landscapes. This analysis, in turn, allows 
us to quantify the degree of efficiency that the network 
has achieved in relation to its operativeness as a global 
transport system. In particular, the major role played by 
interfaces, bridging the communication between the dif- 
ferent percolation components [13 , allowed us to define 
structural efficiency in terms of two complementary as- 
pects: stress or structural load carried by the interfaces 
-which also informs about robustness-, and closeness or 
extent of the direct access to the processing core. We 
use theoretical arguments to propose the conformation 
of maximal structural efficiency and demonstrate by the 
analysis of real networks that biological systems exposed 
to long-term evolutionary pressure may be much closer 
to optimality than information technologies systems at 
an early stage of development dominated by competitive 
forces. 



II. THE ARCHITECTURE OF PERCOLATION 
LANDSCAPES 

Global communication is essential to develop efficient 
collective behavior. In flow networks, represented as di- 
rected complex networks, global connectivity is ensured 
by the presence of architectural elements that allow to 
traverse the network from the input to the output com- 
ponents. These layouts are best rationalized in the frame- 
work of percolation theory, so we call them percola- 
tion landscapes. Characteristic topologies in the perco- 
lated phase denote a global flux that organize in distinct 
linked components comprising macroscopic portions of 
the system (Fig. 1 gives a schematic representation). In 



FIG. 1: Schematic diagram of the main components in the 
percolation landscape of a flow network. The core at the 
center comprises nodes in the SCO and edges within, forming 
the SCE. Nodes in red belong to the IN and the ICE is formed 
by red links. Nodes in blue belong to the OUT and blue 
links form the OCE. Both interfaces, ITF and OTF, appear 
in black. 



the percolated phase, the traditional node percolation 
map [T4l [TSl IT6| recognizes a core structure, the giant 
strongly connected component (SCC), which vertices can 
communicate with each other following directed paths. 
In many real systems this core is a processing unit which 
transfers input to output, and so it is connected to pe- 
ripheral components. The input comes from an affer- 
ent component, the giant in-component (IN), composed 
by all vertices that can reach the SCC but cannot be 
reached from it, and the output goes to an efferent com- 
ponent, the giant out-component (OUT), made of all ver- 
tices that are reachable from the SCC but cannot reach 
it. Secondary structures such as tubes or tendrils could 
also be present [1]. Changing the perspective from nodes 
to edges, this picture is complemented by the edge per- 
colation map [13 , where the number of relevant struc- 
tures increases to five: the edges pure components, ICE, 
OCE, and SCE, that are formed by edges connecting 
nodes within the IN, OUT, and SCC respectively; and 
edges forming the interfaces, ITF and OTF, that bridge 
the peripheral components (IN and OUT, respectively) 
to the core. 

This pattern is obviously further shaped by system de- 
pendent specificities that are the refiection of functional 
demands and evolutionary and/or adaptive forces. In 
particular, the specific conformation of the interfaces de- 
termines the structural efficiency and robustness of the 
network as a global transport system and the potential 
risk of bottleneck effects. 



A. Percolation landscapes of real networks 

We consider here two different information processing 
systems characterized by global transport phenomena: 
one socio-technological, the Internet, in contraposition 
with one biological, the nervous system of the nematode 
worm Caenorhabditis elegans (C. elegans). The node and 
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FIG. 2: Bar diagrams of the values showed in Table 1 detailing 
the percolation landscapes of the ASR and the CEN networks 
as compared with the randomized counterparts. Top charts 
in violet show the node percolation maps and bottom charts 
in orange show the edge percolation maps. 



TABLE I: Statistics for the percolation landscapes of the ASR 
and the CEN networks and their randomized versions (sub- 
script R, values are average ± standard deviation rounded off 
to the first significative figure) . The sizes of the main compo- 
nents are given in absolute number of nodes and links. The 
average degrees are (fc,;) = (ko) — 1.87 and (h) — (ko) = 6.82 
respectively. NP stands for Node Percolation map and EP for 
Edge Percolation map. 

ASR ASRi; CEN CENj; 

IN 20060 18900±800 29 25±2 

NP see 90 880±40 195 219±2 

OUT 17 120±30 55 36±2 

Main 20167 19900±800 279 279±0 

TOTAL 24545 24545 279 279 

ICE 20180 6500±400 9 3±2 

ITF 10833 22000±2000 175 154±8 

SCE 389 2300±200 1322 1490±20 

OTC 226 lOOOilOO 330 230±20 

OCE 12 150±50 36 4±3 

Main 31640 32000±2000 1872 1883±5 

TOTAL 45914 45914 1903 1903 



edge percolation maps of their directed network recon- 
structions are detailed in Fig. 2 and Table 1. 

Their maximally random counterparts are also ana- 
lyzed as null models. In practice, the randomization 
is achieved at the stationary state of a rewiring pro- 
cess that at each time step randomly selects a couple 
of links and exchange their ending points |25j avoiding 
the formation of multiple and self-connections and bidi- 
rectional links while preserving the degree distribution 
P{k). The randomized version would preserve as well 
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degree-degree correlations and higher-order effects which 
correspond to structural constraints ensuring the realiz- 
ability of the network . The comparison of real net- 
works with their randomized counterparts makes there- 
fore possible to determine to which extent the measured 
values are due to global organizing principles and not to 
random assemblages affected by finite-size effects. In this 
work, the randomized counterparts are produced out of 
10 randomized realizations. 



1. Internet customer-provider AS relationships (ASR) 

The Internet is one of the paradigmatic information 
technology and communication networks |17| . From an 
operative point of view, it is composed of thousands of 
Internet Service Providers, usually identified with au- 
tonomous systems (ASs), that operate individual parts 
of the whole infrastructure and engage in contractual re- 
lationships to collectively route traffic through the net- 
work. Such business dependencies jT8] are mappable to a 
directed graph representation of unambiguous customer- 
provider relationships among ASs. 

The directed graph is reconstructed from the 
map 2007-04-02 of inferred AS relationships pro- 



vided by CAID A (http://www.caida.org/data/ active/as-| 
[relationships/). Relationships among ASs are usually 
realized in the form of business agreements, generally 
simplified to customer-provider, peer-to-peer and sibling- 
to-sibling. In a purely directed version of the network, 
where links represent net flow of payments for services 
provided, relations between siblings immediately cancel 
out since they administratively belong to the same orga- 
nization. Peer-to-peer relations are however not trivial 
because they just freely exchange traffic between them- 
selves and their customers but not up in the hierarchy. 
Anyway, we assume here that the later are balanced in 
both directions so as a first approximation we neglect 
them as well. On the other hand, customer-provider re- 
lationships are unambiguously represented by directed 
edges from customer to provider. We are left with a 
purely directed network of 24545 nodes and 45914 di- 
rected links, after removing 4312 (8.55%) peer-to-peer 
and 236 (0.47%) sibhng-to-sibling relations). The in- 
degree distribution is very broad and well described by 
a power law with characteristic exponent 2.1. The out- 
degree distribution is strongly bounded and decays ex- 
tremely fast with a maximum out degree of 24. 

This network presents an extremely asymmetric struc- 
ture at the level of the node percolation map, with a huge 
IN component, a restricted SCC, and an even smaller 
OUT component (Table 1). By comparison, the random- 
ized counterpart is characterized by a similar IN compo- 
nent, but by ten-fold larger (albeit still small in abso- 
lute terms) SCC and OUT components. This informa- 
tion about the node partition should be complemented 
by the analysis of the edge percolation map to provide a 
first glimpse of the different architectural organization of 



the real versus the randomized network. Again the size 
in number of edges of the core and the efferent struc- 
tures (see left graph in Fig. 3 for details about the effer- 
ent components) are qualitatively in accordance with the 
values for the randomized network, despite being smaller. 
However, the organization of the afferent components is 
very different from random. The ICE of the real network 
contains as many edges as nodes in the IN component. 
Moreover, the number of ITF edges connecting the IN 
and SCC components is just half the number of IN nodes: 
on average, thus, there are two IN nodes for every ITF 
edge, which further implies that many nodes in the in- 
component lack direct access to the core. By converse, 
the randomization predicts an ITF double in number of 
edges than actually observed with a correspondingly re- 
duced ICE, so a more shallow IN. 



2. Synaptic neuronal structure of C. elegans ( CEN) 

A different family of information transport systems 
that naturally emerge as archetypical networks are bio- 
logical nervous systems. As for most other complex net- 
works, their structure is intimately related to their func- 
tion and the emergent behavior cannot be understood 
from the mere summation of the individual neuronal ac- 
tions. We focus on the nervous system of the C. elegans 
worm which is practically completely known |19j . 

Network representations of brains display neurons 
as vertices and connection between pairs are present 
whenever a synapse or gap junction has been ob- 
served. We use the updated data set presented in [T5] 
(http://www.wormatlas.org/). The pharyngeal system 
comprises 20 neurons and is almost totally disconnected 
from the rest of the network. It is excluded along un- 
connected neurons, as well as connections of the somatic 
nervous system to non-neural cells. We further restrict to 
chemical synapses excluding gap junctions, very different 
from the previous in nature and function. For simplic- 
ity, polarity or multiplicity of connections are not taken 
into account but directionality is. The synapses are di- 
rected in nature but 233 reciprocal connections has been 
detected (12%). We handled this issue by exploiting the 
imbalance in the number of observed synapses in each 
direction, so that we preserve the directionality of the 
larger number. In this treatment, just 58 of them can- 
cel out (3%). The final set contains 279 nodes and 1903 
links. As reported previously, it turns out to be a small- 
world network [23' with tails of the cumulative distribu- 
tion of degrees for both incoming and outgoing neuronal 
links that have been reported to be well approximated 
by exponential decays fM'. 

Its percolation layout is surprisingly close to random 
organization. In contrast to the Internet, the main struc- 
ture consists of a big core with an OUT twice as large as 
the IN one (see right graph in Fig. 3 for details about the 
afferent components), in accordance to the randomized 
counterpart. The number of edges within the peripheral 



FIG. 3: Left. Efferent structure in ASR. Edges in blue belong to the OCE and connect blue nodes within the OUT. Right. 
Afferent structure in CEN. Edges in red form the ICE and represent connections within red nodes in the IN. In both represen- 
tations, nodes in green belong to the core and edges in black form the corresponding interfaces, OTF and ITF respectively. 



components is extremely small, so that the C. elegans 
nervous system seems to rely on clear input and output 
signals with direct access to the SCC, the computational 
processing core, through well populated interfaces. 



III. INTERFACES AND STRUCTURAL 
EFFICIENCY 

Interfaces play the pivotal role of connecting the IN 
and OUT components to the network core, the SCC com- 
ponent. Setting aside the discussion of wiring costs [20] . 
the efficiency of interfaces at fulfilling such task may be 
loosely defined using complementary measures able to 
capture both the amount of load that the interface edges 
must bear, the "stress", and the extent of the direct ac- 
cess that peripheral component nodes have to the SCC 
through the interface, the "closeness". 



A. Stress as random-walk betweenness 

As elements transported in the system travel the net- 
work, edges are subject to loads that can be characterized 
as betweenness, a topological measure of the number of 
paths between nodes in different components that tra- 
verse those edges. Betweenness is thus a measure of the 
extent to which such edges have control or are under 
stress because of the flow passing through them. Typi- 
cally, betweenness is calculated taking into account only 
shortest-paths between pairs of nodes [21]. Here, we are 
however interested in more realistic situations and as- 
sume that the topological structure is supporting blind 
flow without global knowledge of the system. A more 



appropriate measure is therefore the random-walk be- 
tweenness [22 , that counts all possible routes assuming 
that information wanders at random until it flnds the 
target. Edges with higher random-walk betweenness are 
expected to be more important for the spread of infor- 
mation across the system and, if the load is excessive, 
bottleneck effects could even appear. 

In order to calculate the random-walk betweenness of 
the edges at the interfaces, we slightly modify the origi- 
nal proposal as a centrality measure for vertices The 
percolation landscape is explored by means of two sym- 
metric random walks on the unweighted directed network 
with homogeneous diffusion probabilities and absorbing 
sinks in the nodes of the SCC. Nodes in the IN act as 
sources of diffusive particles -either units of energy, pack- 
ets of information, economic goods, monetary units... - 
which spread from neighbor to neighbor following outgo- 
ing links, each chosen with equal probability among the 
possibilities. The hopping process is stopped whenever 
the diffusive particle arrives to a node in the SCC follow- 
ing a given link in the ITF, which receives the annotation. 
The symmetric process originates particles in the nodes 
of the OUT, which travel backwards following incoming 
links selected with equal probability among the possi- 
bilities, and the diffusion is equally stopped whenever a 
node in the SCC is reached through a particular link in 
the OTF, which receives the annotation. By repeating 
the processes a sufficient number of times for each source 
node it is possible to obtain the probability vector that 
a traveling unit originated at one of the peripheral com- 
ponents uses the edge j in the corresponding interface to 
reach the core. After multiplying by the size of the source 
component in number of nodes Np, the resulting vector 
bjj informs about the structural load that each link in 
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FIG. 4: Upper row panels. Cumulative random-walk be- 
tweenness distribution for the in- and out-interfaces. The 
empirical distributions (symbol lines) are compared against 
those for the randomized null models (dashed lines). Bottom 
row panels. Fraction of nodes remaining in the peripheral 
components of the real networks after removing a fraction p 
of edges of the corresponding interfaces. A targeted removal 
by load in decreasing order (symbol lines) is compared with 
a random deletion (dashed lines). 



interface. At the same time, the results should be again 
validated by investigation of the maximally random coun- 
terpart. In Fig. 4 (upper panels), grey areas denote stress 
regions with loads above 1. Whereas the CEN network 
entirely conforms once more to the randomized predic- 
tion with the practical totality of loads below the thresh- 
old, most edges of the in-interface of the ASR network 
appear to be overstressed, a clear indicator of the vulner- 
ability of the system. The region of loads much below 1 
usually corresponds to peripheral leaf nodes connected to 
multiple core nodes. Apart from a signature of local ro- 
bustness, this diversification could be interpreted as well 
as a quality of being a peripheral spreader or collector of 
flow. 

Finally, the average stress-related structural efficiency 
of an interface can be simply approximated as 



(1) 



that is, the average number of interface edges that medi- 
ate between peripheral nodes and the SCC. This average 
coincides with the inverse of the average betweenness of 
the edges at the interfaces, (-B)j = '^jep^jif-^i- Higher 
values of (fc^) are clearly desirable as peripheral nodes 
would have more routes to the SCC. 



the interface supports. Vector bji/Np corresponds to a 
normalized probability distribution whenever tendrils or 
tubes are not considered. The presence of those appen- 
dices produce cul-de-sac which receive part of the diffu- 
sion unloading partially the interfaces. 

In Fig. 4 (upper panels), we provide the cumula- 
tive distribution Pc{b) = J2b'>b ^i^') values bjj, 
the random-walk betweenness for edge j in interface / 
(I=ITF or OTF), which correspond to the loads of the 
edges at the interfaces of the ASR and CEN networks. 
The cumulated probability density function of the loads 
shows that they are not uniformly distributed for either 
network but have heavy tails denoting large fluctuations, 
with a few links bearing a much higher level of structural 
stress. This heterogeneity is not per se indicating that the 
interface is overstressed. The random-walk betweenness 
is moderately highly correlated with degree [22] mean- 
ing that, in general, vertices with higher degree tend also 
to have higher random- walk betweenness, so that strong 
disorder in the topology could induce spurious hetero- 
geneity in the load distribution. 

In order to assess whether the structural load could 
represent a potential danger of bottleneck formation in 
traffic related processes running on the network, one has 
to define further what is expected as a low load in the 
situation of maximal structural efficiency. We make the 
assumption that such efficiency is reached whenever each 
edge in the interface carries at most a unitary load. This 
gives a simple criteria which makes possible to compare 
different networks but also different links of the same 



B. Stress and robustness 

The loads of the edges at the interfaces are related to 
their robustness, deflned as a measure of the ability of 
the interfaces to communicate different components un- 
der malfunction or failure. In the bottom panels of Fig. 4, 
we show the fraction of nodes remaining connected in the 
peripheral components after the removal of an increasing 
fraction of edges at the corresponding interface. Two dif- 
ferent experiments are performed, the first choosing edges 
according to load in decreasing order and the second se- 
lecting them at random. The results prove that although 
the interfaces seem to be quite robust against random 
failures, the failure of high load edges would disconnect 
a bigger portion of peripheral nodes, thus strongly af- 
fecting the behavior of the system. The CEN and AS 
networks substantially differ in this respect. About 40% 
of interface edges must be removed in CEN before 50% 
of the peripheral nodes are disconnected from the SCC in 
the targeted experiment. By converse, the AS network is 
more delicate because the same degree of disconnection 
is reached by removing just 20% of the interface edges. 

C. Closeness 

The random walk methodology presented above can- 
not discriminate between peripheral conformations with 
different access to the SCC if equal loads are associated to 
interface edges (as a simple example, see tree-like groups 
A and B in Fig.l). The concept of closeness allows us to 
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FIG. 5: Bar diagrams summarizing fine details of tfie peripfi- 
eral components and the interfaces of ASR and CEN networks 
as compared with the randomized null models. Charts in red 
refer to the afferent structure and those in blue to the efferent 
one. The first two bars in each graph refer to the peripheral 
node components detailing the proportions of nodes at dis- 
tance 1 and larger and the proportion of leaf versus non-leaf 
nodes respectively, and the third bar refers to leaf versus non- 
leaf edges at the interfaces. 



TABLE II: Interfaces and peripheral components fine details 
for the ASR and CEN networks and their randomized coun- 
terparts (subscript R, average ± standard deviation rounded 
off to the first significative figure) . Edges at the interfaces and 
nodes at the components are separated into leafs (l) and non- 
leafs (nl), and nodes directly connected to the core {d — 1) 
are distinguished from those at larger distances {d > 1). 
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shed light on the different efficiencies that characterize 
these dissimilar architectures. 

By convention, leaf vertices are those with in-degree 
or out-degree 0, so that they are restricted to belong to a 
peripheral component. In-leaf edges (out-leaf edges) are 
considered as directed links leaving from (pointing to) a 
leaf vertex^?]. Non-leaf edges in the interfaces are the 



TABLE III: Structural efficiency average degrees for ASR and 
CEN and their randomized counterparts (subscript R, values 
are average ± standard deviation rounded off to the first sig- 
nificative figure). Infinite closeness averages come from the 
fact of all peripheral nodes being directly connected to the 
core. 
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ones that ensure the communication from/to nodes not 
directly connected to the core. These non-leaf edges are 
the potentially responsible for bottleneck effects, since 
they service from more than a single IN or OUT node. 
A first estimate of how this topological considerations 
affects efficiency at the structural level is given by the 
closeness average degree, 

{kc) = (2) 

which is the number of interface non-leaf edges available 
for each peripheral node which is not directly connected 
to the sec (thus, with a distance d from the SCC greater 
than I). 

Values for the decomposition of the ASR and CEN 
interfaces and the peripheral components into leafs and 
non-leaf units along with average degree efficiency mea- 
sures as defined in Eq. ([ij and Eq. ^ are shown in Ta- 
ble 2, Table 3, and Fig. 5. In general terms, the higher 
the averages the more structurally efficient the system 
is expected to be. An important imbalance is observed 
between the in and out values for ASR. According to the 
average values, the in-interface presents a certain level 
of inefficiency, with low average degrees combined with 
a low number of leafs, much below random expectations. 
In this situation, potential bottleneck effects are more 
likely. In contrast, the out component shows high lev- 
els of structural efficiency, with the practical totality of 
nodes being root nodes directly connected to the core. 
On the other hand, all peripheral CEN nodes have direct 
access to the core, a signature of high efficiency. 

D. Maximum structural efficiency and the "hairy" 

ball 

Under the requirements of low stress and high close- 
ness, and in the approximation of inexpensive edges, 
maximum efficiency would be realized by a percolation 
landscape structured as a perfect "hairy ball", with all 
the nodes in the peripheral components directly attached 
to the core through leaf edges, each carrying at most a 
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unitary load [25], thus without endangering bottleneck 
effects. Moreover, the interfaces would be robust because 
the failure or malfunctioning of any of its edges would af- 
fect a minimum number of nodes in the peripheral com- 
ponents. Finally, all peripheral nodes would have direct 
access to the core. Any departure from the "hairy ball" 
paradigm would lead to situations in which at least one 
of the two or both requests for structural efficiency, low 
stress and high closeness, are violated to some degree. 

IV. CONCLUSIONS 

Our thorough analysis of percolation landscapes shows 
that the conformation of interfaces plays a central role 
in the performance of complex flow networks as global 
transport systems, governing their efficiency against bot- 
tlenecks and their robustness against failures. We high- 
light that, from the purely structural efficiency perspec- 
tive, a "hairy ball" design would be optimal. Appeal- 
ingly, such behavior may be even displayed by a very close 
to random architecture as seen for the synaptic neuronal 
network of C. elegans. Of the two real systems analyzed 
in this work, this is the one much closer to such opti- 
mality whereas the Internet network presents inefficien- 
cies. These findings point to two, not mutually exclusive, 
interpretations. On the one hand, different adaptation 
dynamics are surely at work: whereas the present struc- 
ture of the C. elegans nervous system tries to optimize 
its collective performance without inter-neuron competi- 
tion, the Internet network emerges, due to its customer- 
provider relations, as a competitive network where it is 



not the global optimization which is sought for but rather 
the individual Internet service provider gain. In this re- 
spect, global efficiency is important only in relation to its 
marketable value. Interestingly, the Internet customer- 
provider network outperforms its randomized version in 
the OTF and OUT components, which describe the ulti- 
mate cash ffow, and underperforms it in the afferent com- 
ponents, where end-users are. On the other hand, evo- 
lution of the worm nervous system might have allowed 
better architectures to emerge, due to its evolutionary 
time-scale (hundred of millions of years) running much 
longer than the time-span of existence of the commercial 
Internet network (slightly more than ten years). 

Clearly, these results only shed light on the basic struc- 
tural ingredients for efficiency and robustness. Indeed, 
several other constraints {e.g. costs of edge deployment 
and maintenance or capacity), are at play which should 
be taken into account for more precise and system spe- 
cific analysis. Yet, percolation landscapes represent a 
first general framework to highlight potential problems 
in a network structure, possibly suggesting specific ac- 
tions to reinforce stressed elements or the redistribution 
of loads so to reduce the risk of bottlenecks and the im- 
pacts of failures. 
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